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Abstract 

We propose a new implementation of a universal set of one- and two-qubit 
gates for quantum computation using the spin states of coupled single-electron 
quantum dots. Desired operations are effected by the gating of the tunneling 
barrier between neighboring dots. Several measures of the gate quality are 
computed within a newly derived spin master equation incorporating decoher- 
ence caused by a prototypical magnetic environment. Dot- array experiments 
which would provide an initial demonstration of the desired non-equilibrium 
spin dynamics are proposed. 
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I. INTRODUCTION 



The work of the last several years has greatly clarified both the theoretical potential and 
the experimental challenges of quantum computation (jT). In a quantum computer the state 
of each bit is permitted to be any quantum-mechanical state of a qubit (two-level quantum 
system). Computation proceeds by a succession of "two-qubit quantum gates" 0, coherent 
interactions involving specific pairs of qubits, by analogy to the realization of ordinary 
digital computation as a succession of boolean logic gates. It is now understood that the 
time evolution of an arbitrary quantum state is intrinsically more powerful computationally 
than the evolution of a digital logic state (the quantum computation can be viewed as a 
coherent superposition of digital computations proceeding in parallel). 

Shor has shown || how this parallelism may be exploited to develop polynomial-time 
quantum algorithms for computational problems, such as prime factoring, which have pre- 
viously been viewed as intractable. This has sparked investigations into the feasibility of 
the actual physical implementation of quantum computation. Achieving the conditions for 
quantum computation is extremely demanding, requiring precision control of Hamiltonian 
operations on well-defined two-level quantum systems, and a very high degree of quantum 
coherence j3j. In ion-trap systems ||, and in cavity quantum electrodynamic experiments 
||, quantum computation at the level of an individual two-qubit gate has been demon- 
strated; however, it is unclear whether such atomic-physics implementations could ever be 
scaled up to do truly large-scale quantum computation, and some have speculated that 
solid-state physics, the scientific mainstay of digital computation, would ultimately provide 
a suitable arena for quantum computation as well. The initial realization of the model that 
we introduce here would correspond to only a modest step towards the realization of quan- 
tum computing, but it would at the same time be a very ambitious advance in the study 
of controlled non-equilibrium spin dynamics of magnetic nanosystems, and could point the 
way towards more extensive studies to explore the large-scale quantum dynamics envisioned 
for a quantum computer. 
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II. QUANTUM-DOT IMPLEMENTATION OF TWO-QUBIT GATES 

In this report we develop a detailed scenario for how quantum computation may be 
achieved in a coupled quantum dot system 0. In our model the qubit is realized as the 
spin of the excess electron on a single-electron quantum dot, see Fig. 1. We introduce 
here a novel mechanism for two-qubit quantum-gate operation that operates by a purely 
electrical gating of the tunneling barrier between neighboring quantum dots, rather than by 
spectroscopic manipulation as in other models. Controlled gating of the tunneling barrier 
between neighboring single-electron quantum dots in patterned two-dimensional electron- 
gas structures has already been achieved experimentally using a split-gate technique ||. If 
the barrier potential is "high", tunneling is forbidden between dots, and the qubit states 
are held stably without evolution in time (t). If the barrier is pulsed to a "low" voltage, the 
usual physics of the Hubbard model @ says that the spins will be subject to a transient 
Heisenberg coupling, 

H s (t) = J(t)S 1 -S 2 , (1) 

where J(t) = —^p- is the time-dependent exchange constant |L0] which is produced by the 
turning on and off of the tunneling matrix element t (t). Here u is the charging energy of a 
single dot, and Si is the spin-1/2 operator for dot i. 

Eq. (P will provide a good description of the quantum-dot system if several conditions 
are met: 1) Higher-lying single particle states of the dots can be ignored; this requires 
AE ^> kT, where AE is the level spacing and T is the temperature. 2) The time scale t s for 
pulsing the gate potential "low" should be longer than h/ AE, in order to prevent transitions 
to higher orbital levels. 3) u > to(t) for all t; this is required for the Heisenberg-exchange 
approximation to be accurate. 4) The decoherence time T _1 should be much longer than the 
switching time r s . Much of the remainder of the paper will be devoted to a detailed analysis 
of the effect of a decohering environment. We expect that the spin-1/2 degrees of freedom 
in quantum dots should generically have longer decoherence times than charge degrees of 
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freedom, since they are insensitive to any environmental fluctuations of the electric potential. 
However, while charge transport in such coupled quantum dots has received much recent 



attention |Tl|,|8| we are not aware of investigations on their non-equilibrium spin dynamics 
as envisaged here. Thus, we will carefully consider the effect of magnetic coupling to the 
environment. 

If T _1 is long, then the ideal of quantum computing may be achieved, wherein the 
effect of the pulsed Hamiltonian is to apply a particular unitary time evolution operator 
U 3 (t) = Texp{-i Jq H s (t')dt'} to the initial state of the two spins: |#(t)) = U s \^(0)). The 
pulsed Heisenberg coupling leads to a special form for U s : For a specific duration t s of the 
spin-spin coupling such that / dtJ{t) = J t s = it (mod 2n) fL2f , U s (J t s = n) = U sw is the 
"swap" operator: if labels the basis states of two spins in the S^-basis with i,j = 0, 1, 
then U sw \ij) = Because it conserves the total angular momentum of the system, U sw 
is not by itself sufficient to perform useful quantum computations, but if the interaction is 
pulsed on for just half the duration, the resulting "square root of swap" is very useful as 
a fundamental quantum gate: for instance, a quantum XOR gate is obtained by a simple 
sequence of operations: 

Uxor = e^e'^uL^uL , (2) 

where e mS l etc. are single-qubit operations only, which can be realized e.g. by applying 
local magnetic fields (see Sec. 1 [ [B|) ||13|| . It has been established that XOR along with 



single-qubit operations may be assembled to do any quantum computation 0. Note that 
the XOR of Eq. (fj) is given in the basis where it has the form of a conditional phase-shift 
operation; the standard XOR is obtained by a simple basis change for qubit 2 . 

III. THE MASTER EQUATION 

We will now consider in detail the non-ideal action of the swap operation when the two 
spins are coupled to a magnetic environment. A new master equation model is obtained that 



4 



explicitly accounts for the action of the environment during switching — to our knowledge, 
the first treatment of this effect. We use a Caldeira-Leggett type model in which a set 
of harmonic oscillators are coupled linearly to the system spins by H int = XJ2i=i2^i ' 
Here, b{ = J2 a 9a( a a,ij + a a ,ij) * s a fluctuating quantum field whose free motion is governed 
by the harmonic oscillator Hamiltonian H B = ^oo^a^^a^ij, where a a> iJa a ,ij are bosonic 
creation/annihilation operators (with j = x, y, z), and to 7 * are the corresponding frequencies 
with spectral distribution function Jij(u) = ^ ) 2 <5(w — oj a ) [fL4| . The system and 

environment are initially uncorrelated with the latter in thermal equilibrium described by 
the canonical density matrix p B with temperature T. We assume for simplicity that the 
environment acts isotropically and is equal and independent on both dots. We do not 
consider this to be a microscopically accurate model for these as-yet-unconstructed quantum 
dot systems, but rather as a generic phenomenological description of the environment of a 
spin, which will permit us to explore the complete time dependence of the gate action on 
the single coupling constant A and the controlled parameters of H s (t) fl"5|| . 



A. Swap gate 

The quantity of interest is the system density matrix p(t) = TrBp{t) which we obtain by 
tracing out the environment degrees of freedom. The full density matrix p itself obeys the 
von Neumann equation, 

p{t) = -i[H,p] = -iCp, (3) 

where 

£ = £ s (t) + £ int + £ B (4) 
denotes the Liouvillian |TB[ corresponding to the full Hamiltonian 

H = H s (t)+H mt + H B . (5) 
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Our goal is to find the linear map (superoperator) V(t) which connects the input state of the 
gate p = pit = 0) with the output state pit) after time t > r s has elapsed, p(t) = V(t) p . 
V(t) must satisfy three physical conditions: 1) trace preservation Tr s Vp = 1, where Tr s 
denotes the system trace, 2) Hermiticity preservation (VpY = Vp, and 3) complete positivity, 
(V®1b)p > 0. Using the Zwanzig master equation approach |TB[ we sketch the derivation for 
V in the Born and Markov approximation which respects these three conditions. The novel 
situation we analyze here is unusual in that H s is explicitly time-dependent and changes 
abruptly in time. It is this fact that requires a separate treatment for times t < t s and 
t > r s . To implement this time scale separation and to preserve positivity it is best to start 
from the exact master equation in pure integral form, 

p(t) = U s (t, 0)p - fda [ a dTU s {t, a)M(a, r)p(r), (6) 
Jo Jo 

where 

Ui(t,t') =Tex P {-i f drCiir)}, (7) 
Jt' 

where i — s, B, int, or q. Here q indicates the projected Liouvillian 

C q = (l-P)C = (l-p B Tr B )C. (8) 

Also, the "memory kernel" is 

M{a, t) = Tr B C int l4 q (a, T)C int p B . (9) 

We solve (^|) in the Born approximation and for t ^> t s . To this end the time integrals 
are split up into three parts, 1) < r < a < r s < t, 2) < r < r s < o < t, and 3) 
0<r s <T<a<t. Keeping only leading terms in r s we retain the contribution from 
interval 2) as it is proportional to t s , whereas we can drop interval 1) which leads to higher- 



order terms. But note that terms containing J r s must be kept to all orders |Hq| . Interval 
3) is independent of r s . 

Rewriting the expressions and performing a Born approximation (i.e., keeping only lowest 
order terms in A 2 ) with subsequent Markov approximation we find, for t > r s , 



V(t) = e-^- T ^U s {T s ) (1 - /C 2 ) , (10) 
where Ui{r s ) = Ui(r s , 0), /C2 describes the effect of the environment during the switching, 



£ a = Wfc) / dr / dt 



r s /-co 

dr / 

jO 

x Tr B C in tU s {T)U B {t)C int p B U s {T s - r), (11) 



while 
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/C 3 = / dt Tr B C int U B (t)C int p B (12) 

JO 

is independent of if s . We also note that £4(1 — /C2) has a simple interpretation as being the 
"transient contribution" which changes the initial value p$ ai t = § to U s (t s )(1 — A^po a ^ 
t = t s . We show in the Appendix that, to leading order, our superoperator V indeed satisfies 
all three conditions stated above, in particular complete positivity. Such a proof for spins 
with an explicit time-dependent and direct interaction, Eq. ([[]), is not simply related to the 
case of a master equation for non-interacting spins (and without explicit time-dependence) 
considered in the literature (see for example P~T|JT5| ] ) . We also note that the above Born 
and Markov approximations could also be introduced in the master equation in the more 
usual differential-integral representation. However, it is well-known from studies in non- 
interacting spin problems [18 that in this case the resulting propagator is in general no 
longer completely positive. 

Next, we evaluate the above superoperators more explicitly, obtaining 

K 2 p = (r + iA) ]T r dr[Si(T s ), -SiWp) + h.c. , (13) 

i Jo 

JC 3 p = T(3p-2j2S l p-S l ) , (14) 



where T, A are real and given by 

\ 2 /-oo roc 

T = — / dt duJiuo) cosiut) coth(— — — ) , (15) 
7r Jo jo 2k B T 

/ dt dujJ(uj) sin(cjt) . (16) 
Jo Jo 



A = 

7T 



In our model, the transverse and longitudinal relaxation or decoherence rates of the system 
spins are the same and given by T. For instance, for ohmic damping with J{uj) = rjcu, 
we get T = \ 2 i]kBT and A = \ 2 rjuj c /Ti, with u c some high frequency cut-off. Requiring 
for consistency that Tr s , At s <C 1 we find that /C 2 is in fact a small correction. However, 
we emphasize again that, to our knowledge, this is the first time that any analysis of this 
/C 2 term, describing the action of the environment during the finite time that the system 
Hamiltonian is switched on, has been given. 

For further evaluation of V we adopt a matrix representation, defined by 

Vab\cd = (e ab , Ve cd ) = Tr s e\ b Ve cd , (17) 

where {e a b|a, b = 1, .., 4} is an orthonormal basis, i.e. (e a b, e cd ) = S ac 5b d . In this notation we 
then have 

P(t)ab = Vab\cd(po)cd (18) 

c,d 

with V being a 16 x 16 matrix. 

Note that /C 2 ,3 and U s are not simultaneously diagonal. However, since K,s(l,Si) = 
2r(0, Si) we see that exp{ — (t — r s )/C 3 } is diagonal in the "polarization basis" {e p ab = 
e a e h e \,..,4 = (l/y/2,V2S?,y/2S?,\/2S?),i = 1,2}, while C s and thus U s are diagonal 
in the "multiplet basis" {e^ = \a)((3\,a,(3 = 1,...,4; |1) = (|01) - |10))/v^,|2) = 
(|01> + \10))/V2, |3) = |00), |4) = |11)}, with 

u s {tu wp , = LV^"^ 1 , (19) 

where E™ = — 3J /4, = J /A are the singlet and triplet eigenvalues. Lastly, /C 2 

is most easily evaluated also in the multiplet basis; after some calculation we find that 
/C 2 = K, d 2 - ICf, with 

(fci)a(3\j5 = Y^^-y( S \^\ a ') ' ( a '\Si\P)K>a>\5(3 + ^ S (a\Si\a'} ■ (a'\Si |7>^a'a'| 7 a] , (20) 
i,a' 

(JCfU b s = I>$l7> • (8\$i\P) (kapw + (k MSl )*) . (21) 

i 
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Here 

k aphs = (r + iA)e l{E ™- E ™ )T ° ^ dre^- E ^ 

= 2^[ Fc <5/3 ~ As <5/3 + i ( Ts Sf3 + Ac 5/3 )][s a7 + l(l - C Q7 )], 

Cij = cos^TsUij), s i:j = sinfoafy), LUij = E™ - Ef. 
Using the above matrix notation, we can write explicitly 

V«f3hs= E (^| a/3 )*(e- ( *- Ts) ^) a%fe C , afc | a ^,e- irs( ^-^ ) (l-/C 2 ) a ^ l75 , (22) 

a,b,a',p' 

where Cab\ a p = (e^ 6 , e^) is the unitary basis change between the polarization and the 
multiplet basis. 

B. One-bit gates 

We now repeat the preceding analysis for single-qubit rotations such as e l % s * as required 
in Eq. (0). Such rotations can be achieved if a magnetic field Hi could be pulsed exclusively 
onto spin i, perhaps by a scanning-probe tip. An alternative way, which would become 
attractive if further advances are made in the synthesis of nanostructures in magnetic semi- 



conductors is to use, as indicated in Fig. 1(a), an auxiliary dot (FM) made of an 
insulating, ferromagnetically-ordered material that can be connected to dot 1 (or dot 2) by 
the same kind of electrical gating as discussed above 0. If the the barrier between dot 1 
and dot FM were lowered so that the electron's wavefunction overlaps with the magnetized 
region for a fixed time t s , the Hamiltonian for the qubit on dot 1 will contain a Zeeman term 
during that time. For all earlier and later times the magnetic field seen by the qubit should 
be zero; any stray magnetic field from the dot FM at neighboring dots 1, 2, etc. could be 
made small by making FM part of a closure domain or closed magnetic circuit. 
In either case, the spin is rotated and the corresponding Hamiltonian is given by 

r'dtHf = j2^sSf, (23) 
Jo i=i 
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with uji = g/^BH?, where we assume that the H-field acting on spin i is along the z-axis. The 
calculation proceeds along the same line as the one described above: Just as in Eq. (|IUD , 
the expression obtained for the superoperator is 

V„{t) = e-^**Uf (t s ) (1 - K* ) . (24) 

/C3 is exactly the same as before, Eq. (0). Uf(r s ) is again given by Eq. (|7|) with the mod- 
ification that the Liouvillian (see Eq. (^)) corresponding to the magnetic field Hamiltonian 
of Eq. ( |23| ) is used rather than that for the exchange Hamiltonian H s (Eqs. (0) and ([!])). 
The explicit matrix representation is 

(U? (r a ))„| rV = 6 rr ,5 ss ,e-^l^- El >°. (25) 

Here we are employing another basis, the S^-basis for the two spins {e z rs = \r)(s\, r,s — 
1,2,3,4; \s) = |00), |01), |10), |11)}. The energies are 

K 2 } = Kw) = -1,1,-1}- (26) 

The K-2 calculation also proceeds as before (see Eq. using the new Hamiltonian; the 

result is K,2 — /C 2 ' — /C 2 , with 

(/Cf' d ) r ^ u = £ (Mu$|r') • (r'|^|s)(^ V |„ s )* + S su {r\Sy) • (r'|^|t)^ v|tr ) , (27) 

i,r' 

{Kf^Utu = £<r|S^> • <«$l a > + J*) • (28) 

i 

Here 

fc* a|tu = (r + zA)e*( B «- £, i) r * / rs dTe^-^ 

= 2^[ r C - + *( r C + A C)]fe + i(i - 4*)L ( 29 ) 

The i? fc, s are from Eq. (^). Finally, the explicit matrix form for V H may be written 

V%>\a'V = £ (e- (t - rs)/C3 )a% 6 (^r S |a 6 )*e-^'=i^^-^(l - /Cf )„| rV A-V|a'V , (30) 

r,s,r',s' 

where D rs \ ab = (e* 8 , e^ fe ) is now the unitary basis change between the basis and the 
polarization basis. 
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C. Numerical study for swap gate and XOR gate 



Having diagonalized the problem, we can now calculate any system observable; the re- 
quired matrix calculations are involved and complete evaluation is done with Mathematica. 
We will consider three parameters (s, F, and P in Fig. 2) relevant for characterizing the gate 
operation. We first perform this analysis for the swap operation introduced above. 

The swap operation would provide a useful experimental test for the gate functionality: 
Let us assume that at t — spin 2 is (nearly) polarized, say, along the z-axis, while spin 
1 is (nearly) unpolarized, i.e. po — (1 + 2S'f)/4. This can be achieved, e.g., by selective 
optical excitation, or by an applied magnetic field with a strong spatial gradient. Next we 
apply a swap operation by pulsing the exchange coupling such that Jqt s = it, and observe 
the resulting polarization of spin 1 described by 

(sm) = |v(t) 41 , 14 , (si) 

where V is evaluated in the polarization basis. After time t s spin 1 is almost fully polarized 
(whereas spin 2 is now unpolarized) and, due to the environment, decays exponentially with 
rate of order T. To make the signal Eq. ([31]) easily measurable by conventional magnetom- 
etry, we can envisage a set-up consisting of a large array of identical, non-interacting pairs 
of dots as indicated in Fig. 1(b). 

To further characterize the gate performance we follow Ref. |20| and calculate the gate 
fidelity F = (^o\W(r s )p(t)\xjjo), and the gate purity P = Tr s [p(t)] 2 , where the overbar means 
average over all initial system states |^ )- Expressing V in the multiplet basis and using 
trace and Hermiticity preservation we find 

F(t) = \ + ^ReE V aa \ aa +EV^^], (32) 
p (t) = ^ E l\Vkk'\u\ 2 + £(V fefe '|nV fc V to . + |V fefc '|d 2 )] (33) 

i,k,k' j 

(in fact, the expression for P(t) holds in any basis). Evaluations of these functions for 
specific parameter values are shown in Fig. 2. For the parameters shown, the effect of the 
environment during the switching, i.e. /C2 in Eq. (|10|), is on the order of a few percent. 
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The dimensionless parameters used here would, for example, correspond to the following 
actual physical parameters: if an exchange constant Jq = 80/ieV ~ IK were achievable, then 
pulse durations of t s « 25ps and decoherence times of T" 1 « 1.4ns would be needed; such 
parameters, and perhaps much better, are apparently achievable in solid-state spin systems 

As a final application, we calculate the full XOR by applying the corresponding super- 
operators in the sequence associated with the one on the r.h.s. of Eq. @. We use the same 
dimensionless parameters as above, and as before we then calculate the gate fidelity and 
the gate purity Some representative results of this calculation are plotted in the inset of 
Fig. 2(b). To attain the 7r/2 single-bit rotations of Eq. (|2|) in a t s of 25ps would require a 
magnetic field H « 0.6T, which would be readily available in the solid state. 

IV. DISCUSSION 

As a final remark about the decoherence problem, we would note that the parameters 
which we have chosen in the presentation of our numerical work, which we consider to be 
realistic for known nanoscale semiconductor materials, of course fall far, far lower than the 
0.99999 levels which are presently considered desirable by quantum-computation theorists 
H[; still, the achievement of even much lesser quality quantum gate operation would be 
a tremendous advance in the controlled, non-equilibrium time-evolution of solid-state spin 
systems, and could point the way to the devices which could ultimately be used in a quantum 
computer. 

Considering the situation more broadly, we are quite aware that our proposal for quan- 
tum dot quantum computation relies on simultaneous further advances in the experimen- 
tal techniques of semiconductor nano-fabrication, magnetic semiconductor synthesis, single 
electronics, and perhaps in scanning-probe techniques. Still, we also strongly feel that such 
proposals should be developed seriously, and taken seriously, in the present, since we believe 
that many aspects of the present proposal are testable in the not-too-distant future. This 
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is particularly so for the demonstration of the swap action on an array of dot pairs. Such a 
demonstration would be of clear interest not only for quantum computation, but would also 
represent a new technique for exploring the non-equilibrium dynamics of spins in quantum 
dots. 

To make the quantum-dot idea a complete proposal for quantum computation, we need 
to touch on several other important features of quantum-computer operation. As our guide- 
line we follow the five requirements laid out by one of us [Q: 1) identification of well-defined 
qubits; 2) reliable state preparation; 3) low decoherence; 4) accurate quantum gate oper- 
ations; 5) strong quantum measurements. Items 1, 3, and 4 have been very thoroughly 
considered above. We would now like to propose several possible means by which require- 
ments 2 and 5, for state preparation (read in) and quantum measurement (read out), may 
be satisfied. 

One scheme for qubit measurement which we suggest involves a switchable tunneling (T 
in Fig. 1(a)) into a supercooled paramagnetic dot (PM). When the measurement is to be 
performed, the electron tunnels (this will be real tunneling, not the virtual tunneling used for 
the swap gate above) into PM, nucleating from the metastable phase a ferromagnetic domain 
whose magnetization direction could be measured by conventional means. The orientation 
(8, <p) of this magnetization vector is a "pointer" which measures the spin direction; it is a 
generalized measurement in which the measurement outcomes form a continuous set rather 
than having two discrete values. Such a case is covered by the general formalism of positive- 



operator- valued (POV) measurements [21]. If there is no magnetic anisotropy in dot PM, 
then symmetry dictates that the positive measurement operators would be projectors into 
the over-complete set of spin- 1/2 coherent states 

\0,<p) = cos ^|0) sin ^|1). (34) 

_ _ 

A 75%-reliable measurement of spin-up and spin-down is obtained if the magnetization 
direction (8, 0) in the upper hemisphere is interpreted as up, and in the lower hemisphere 
as down; this is so simply because 
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^/^I<0|M>| 2 = ~- (35) 

Here "£/" denotes integration over the upper hemisphere and 2ir is the normalization constant 
for the coherent states. 

Another approach which would potentially give a 100% reliable measurement requires a 
spin-dependent, switchable "spin valve" tunnel barrier (SV) of the type mentioned e.g. in 
Ref. ||22|| . When the measurement is to be performed, SV is switched so that only an up-spin 



electron passes into semiconductor dot 3. Then the presence of an electron on 3, measured 
by electrometer £, would provide a measurement that the spin had been "up." It is well 
known now how to create nanoscale single-electron electrometers with exquisite sensitivity 
(down to 10~ 8 of one electron) f23|j . 

We need only discuss the state-preparation problem briefly. For many applications in 
quantum computing, only a simple initial state, such as all spins up, needs to be created. 
Obviously, such a state is achieved if the system is cooled sufficiently in a uniform ap- 
plied magnetic field; acceptable spin polarizations of electron spins are readily achievable 
at cryogenic temperatures. If a specific arrangement of up and down spins were needed 
as the starting state, these could be created by a suitable application of the reverse of the 
spin-valve measurement apparatus. 
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APPENDIX A: COMPLETE POSITIVITY OF TIME-EVOLUTION 

SUPEROPERATOR V 



Here we sketch the proof that the superoperator V in Eq. (10) is completely positive. 
We analyze the /C3 term first. We write 

N 

Af^oo V N 



e~*= Jim fl-^3) . (Al) 



It is sufficient to prove that the infinitesimal operator is completely positive. It is straight- 
forward to show, using Eq. (|14]), that 

1 - £fc 3 ) P = Zl-pZ 3 + ((r/N) 2 ) . (A2) 
Here Z 3 is the seven-component vector operator 



2 -^M^ B J- < A3) 



where 

B = (B 1 ,...,B 6 ) = V2T (S U S 2 ). (A4) 

Note that for this case B\ = Bk and Y%=i ^\Bk = 3I\ 

We recall that it is easy to prove that any superoperator S of the form 

Sp = Z^pZ (A5) 

as in the first term of Eq. ( \K*%) is completely positive. Indeed, considering its action on any 
state vector of the system plus environment 0, and taking a positive p we get 

(<f>,Sfxf>) = ((f), Z ] pZcj)) = (Z<f>,pZ(j>) > V0. (A6) 

Next we consider the I — KL2 term of Eq. (|l^). Starting from Eq. (|13D , we put this term 
in a form corresponding to the completely-positive form Eq. (|A5|). We find 

(1 - K 2 )p = Z\ - P Z 2 + (A 4 , rl (A 2 r s ) 2 ) , (A7) 
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with Z 2 being the vector operator 

Z 2 = (l + Y ] ■ X\X -Y x ) (A8) 

with 

X = -{T + iA)(S 1 (r s ),S 2 {r s )), (A9) 
Y= f Ts driSxir),^)). (A10) 







So, from the same arguments as above, Eq. (|A7|) establishes that 1 — /C 2 is completely 
positive up to the order of accuracy discussed in the text. 

Finally, we note that the other two general conditions for a physical superoperator also 
follow immediately: Trace preservation of V follows from the fact that a Liouvillian £ appears 
to the left in the basic equations for /C 2 , Eq- (PD, and /C3, Eq. (12) . Trace preservation is 



also reflected in the fact that Z 2 ■ Z 2 = 1 and Z s ■ Z\ = 1 to leading order. The form Eq. 
( |A5| ) also obviously preserves Hermiticity of the density operator; this is also clear from the 
forms of Eqs. (0) and (0). 
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FIGURES 




FIG. 1. a) Schematic top view of two coupled quantum dots labeled 1 and 2, each contain- 
ing one single excess electron (e) with spin 1/2. The tunnel barrier between the dots can be 
raised or lowered by setting a gate voltage "high" (solid equipotential contour) or "low" (dashed 
equipotential contour). In the low state virtual tunneling (dotted line) produces a time-dependent 
Heisenberg exchange J(t). Hopping to an auxiliary ferromagnetic dot (FM) provides one method 
of performing single-qubit operations. Tunneling (T) to the paramagnetic dot (PM) can be used 
as a POV read out with 75% reliability; spin-dependent tunneling (through "spin valve" SV) into 
dot 3 can lead to spin measurement via an electrometer £. b) Proposed experimental setup for 
initial test of swap-gate operation in an array of many non-interacting quantum-dot pairs. Left 
column of dots is initially unpolarized while right one is polarized; this state can be reversed by a 
swap operation (see Eq. fl3l"D). 
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FIG. 2. a) Swap polarization s = 2(Sf(t)} (see Eq. (|3l])), gate fidelity F, and gate purity P 
vs. Tt for "swap" using parameters Jqt s = tt, Tt s = 0.017, and Ar s = —0.0145. b) The same 
for XOR obtained using the four operations in Eq. (|2|) (the final two single-spin operations done 
simultaneously). The same parameters and scales as in (a) are used; the pulse-to-pulse time is 
taken to be 3t s . Tt is measured from the end of the fourth pulse. 
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